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Abstract. We introduce a stochastic halo formation model 
to compute the early chemical enrichment of the interstellar 
medium (ISM) of the halo. Contrary to 1-zone chemical evo- 
lution models, we are able to resolve local inhomogeneities in 
the ISM caused by single core-collapse supernovae. These in- 
homogeneities lead to different element abundance patterns in 
very metal-poor stars, which can be seen as scatter in the abun- 
dances of halo stars with metallicities [Fe/H] < 2.0. 

The early chemical evolution of the halo proceeds in differ- 
ent enrichment phases: At [Fe/H] < —3.0, the halo ISM is un- 
mixed and dominated by local inhomogeneities caused by in- 
dividual core-collapse supernova (SN) events. For metallicities 
[Fe/H] > — 2.0 the halo ISM is well mixed, showing an element 
abundance pattern integrated over the initial mass function. In 
the range —3.0 < [Fe/H] < —2.0 a continuous transition from 
the unmixed to the well mixed ISM occurs. 

For some elements (Si, Ca, Eu), the scatter in the element- 
to-iron ratio [El/Fe] of metal-poor halo stars can be reproduced. 
Stellar yields of other elements predict a scatter which, com- 
pared to the observations, is too large (O, Mg) or too small 
(Ni). Cr and Mn show a decreasing trend for lower metallic- 
ities, which can not be explained by metallicity independent 
yields, provided that the mixing of the ejecta with the interstel- 
lar medium does not depend on progenitor mass. This demon- 
strates the need for revised, self-consistent SN yields. 

Finally, we discuss the metallicity distribution in the model. 
Compared to the 28 very metal-poor stars observed with metal- 
licities in the range —4.0 < [Fe/H] < —3.0, no star is known 
with confirmed metallicity [Fe/H] < —4.0, while our model 
predicts 5 ± 2 stars with [Fe/H] < —4.0. These should be 
present if the halo ISM started at primordial metallicities and 
no pre-enrichment by population III stars occurred. 

Key words: ISM: abundances - Physical processes: Nucle- 
osynthesis - Galaxy: abundances - Galaxy: halo 



1. Introduction 

The low metal abundances and high peculiar velocities of halo 
stars indicate that the halo is an old, if not the oldest component 
of the Milky Way. Age determinations of globular clusters and 



halo field stars point to an age of 14-15 billion years with no de- 
tectable age gradient with galactocentric distance (Harris et al. 



1997). The halo has therefore special significance for the for- 



mation of the Milky Way. There are two quantities which play 
important roles in the investigation of the formation of the halo: 
the orbits of halo stars and their chemical composition. Since 
halo stars form a collisionless system, their orbits contain in- 
formation about the dynamics at the time of sta r form ation and 
thus th e form ation of the halo (e.g. Carney et al. 1996 ; Chiba & 
Yoshii 1998 ). Information about the chemical composition of 
the interstellar medium (ISM) of the halo is more direct. The 
element abundances observed in low mass halo stars directly 
reflect the chemical abundances and the chemical inhomogene- 
ity of the ISM during halo formation (Mc William |1997[ ). 

Examinations of element abundance ratios as function of 
metallicity [Fe/H] show that star- to-sta r differences rise with 
decreasing metallicity (Ryan et al. 1996 ). Most of the chemical 
elements are ejected during supernovae Type II (SNe II) explo- 
sions. The enrichment of the halo depends on how many SNe II 
explode and how effectively the ejected gas is mixed with the 
surrounding ISM. If the ejected metals are distributed over a 
large volume, a spatially homogeneous enrichment takes place. 
If the mixing volume is small, the ISM in the vicinity of a core- 
collapse supernova (SN II) is highly enriched, while large parts 
of the halo gas remain metal-poor. In this case the ISM is chem- 
ically highly inhomogeneous and newly formed stars are of dif- 
ferent chemical composition, depending on where they form. In 
this scenario one should moreover expect that the metal-poorest 
stars have chemical compositions corresponding to the stellar 
yields of single SNe II (Ryan et al. 



1996) 



In this paper we present a stochastic chemical evolution 
model and investigate the inhomogeneous enrichment of the 
halo ISM. The description of the model is given in Sect. ^, fol- 
lowed by an overview of theoretical SN II yields and their un- 
certainties in Sect. J|. The employed observational data is pre- 
sented in Sect. |j. The results of our model and the conclusions 
are given in Sect. ^ and [| respectively. 

2. The Model 

Very metal-poor halo stars show a great diversity in their el- 
ement abundances and therefore a scatter in their element- 
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to-iron ratios [El/Fe] of order 1 dex. This scatter gradually 
decreases at higher metallicities until a mean element abun- 
dance is reached which corresponds to the [El/Fe] ratio of the 
stellar yields integrated over the initial mass function (IMF). 
The aim of our stochastic halo formation model is to under- 
stand the trends seen in the observations and to investigate how 
the metal-poor interstellar medium (ISM) in the halo evolves 
chemically. 

Our fully 3D-code, contrary to 1-zone chemical evolution 
models, enables us to resolve local inhomogeneities in the ISM 
with a spatial resolution of 50 pc. All in all, we model a volume 
of (2.5 kpc) 3 , divided into 50 3 cells. Every cell of our grid con- 
tains detailed information about the enclosed ISM and the mass 
distribution of stars. We consider simultaneously the evolution 
of nine elements, the a-elements O, Mg, Si and Ca, the iron- 
peak elements Cr, Mn, Fe and Ni and the r-process element Eu. 

Our initial conditions assume a halo ISM consisting of a 
homogeneously distributed single gas phase with primordial 
abundances and a density of 0.25 particles per cm 3 , which 
gives a total mass of about 10 8 Mq in a volume of (2.5 kpc) 3 . 
We adopt a constant time-step of 10 6 years since it has to be 
longer than the dynamical evolution of a supernova (SN) rem- 
nant and shorter than the lifetime of the most massive stars. At 
each time-step 20 000 cells are chosen randomly and indepen- 
dently of each other and of the state of the enclosed ISM. Each 
selected cell may create a star with a probab ility p roportional 
to the square of the local ISM density (Larson 1988 ). The num- 
ber of stars formed per time-step is the product of the number 
of cells tested with the probability of star formation in each 
cell. Various combinations of these parameters are possible to 
achieve a given SFR; the choice of 20 000 cells proved compu- 
tationally convenient. The absolute value of the SFR influences 
the time-scale of the enrichment process (cf. Sect. 5.1 ), but not 
the evolution of [El/Fe] -ratios as function of [Fe/H]. Therefore, 
the main results of this paper are insensitive to the values of 
these parameters. 

The mass of a newly formed star is chosen randomly from 
a Salpeter IMF. The lower and upper mass limits of the IMF are 
taken to be 0.1 Mq and 50 Mq, respectively. About 5000 stars 
are formed on average during each step. Newly born stars in- 
herit the abundance pattern of the ISM out of which they form, 
carrying therefore information about the state of the ISM at the 
place and time of their birth. To determine the lifetime of a star 
an approximation to the metallicity dependent mass-lifetime 
relation of the Geneva Stellar Evolution and Nucleosynthesis 



the neighbouring ISM. Stellar yields are taken from Thiele- 
mann et al. ( 1996j ) and Nomoto et al. (1997) for all elements 
except Eu. Since there are no theoretical predictions of stellar 
Eu yields, we use the indirectly deduced yields of Tsujimoto & 



Group (cf. Schaller et al. 1992; Schaerer et al. 1993a; Schaerer 
et al. 1993b| ; Charbonnel et al. 1993) is used, given by 



log (T) = (3.79 + 0.24 • Z) - (3.10 + 0.35 ■ Z) log (M) 
+ (0.74 + 0.11 ■ Z) log 2 (M), 

where T is the lifetime in units of 10 6 yr, Z the metallicity in 
units of solar metallicity Zq and M the mass in units of solar 

masses Mq. 

Stars in a range of 10 — 50 Mq will explode as core- 
collapse supernovae (SNe II), resulting in an enrichment of 



Shigeyama (1998) which assume that r-process elements origi- 
nate from SNe II (see the discussion of stellar yields in Sect. |]). 
We linearly interpolate the stellar yields given in these papers, 
since we use a finer mass-grid in our simulation. For SNe with 
masses below 13 Mq stellar yields are not available. Since the 
nucleosynthesis models show declining yields towards low pro- 
genitor masses, we have for the interpolation arbitrarily set the 
yields of a 10 Mq SN to one thousandth of those of a 13 Mq 
SN. The interpolation gives IMF averaged values of the [El/Fe] 
ratios, which are in good agreement with the observed mean 
values of metal-poor stars in all elements except Ca, which 
shows a [Ca/Fe] ratio that is about 0.3 dex lower than the ob- 
served mean. We do not include supernovae (SNe) of Type la, 
since we are only interested in the very early enrichment of the 
halo ISM, which is dominated by SNe of Type II. 

Intermediate mass stars will evolve to planetary nebulae, 
returning only slightly enriched material in the course of their 
evolution. This locally influences the enrichment pattern of 
the gas, since metal-poor material is returned into the evolved 
and enriched ISM. It will not change the element abundances 
[El/H] significantly, but can affect the local element-to-iron ra- 
tios [El/Fe] considerably. Low mass stars do not evolve signif- 
icantly during the considered time. In our model, they serve to 
lock up part of the gas mass, affecting therefore the local ele- 
ment abundances [El/H] in the ISM. 

Since the explosion energy of a core-collapse supernova 
(SN II) depends only slightly on the mass of its progenitor 
(Woosley & Weaver|1995|; Thielemann et al. |1996[), every SN II 



sweeps up a constant mass of about 5 x 10 4 Mq of gas (Ryan et 
al. 1996; Shigeyama & Tsujimoto 1998). In our model the ra- 



dius of the SN remnant then is computed from the local density 
of the ISM and lies typically between 100 pc and 200 pc. The 
ejecta of the SN II and all the swept up, enriched material are 
condensed in a spherical shell which is assumed to be chemi- 
cally well mixed. The material in the shell subsequently mixes 
with the ISM of the cells where the expansion of the remnant 
stopped. The interior of the remnant, where all the material was 
swept up, is assumed to be filled with about 5 Mq of dilute gas 
from the SN event with the corresponding metal abundances. 
This gas is unable to form stars until it is swept up by another 
SN event and mixed with the surrounding ISM. Thus this ma- 
terial contributes to the enrichment only after some delay. 

The star formation rate of cells influenced by the remnant 
will rise, since their density is higher than the average density 
of a cell and the probability to form a star is assumed to be 
proportional to the square of the local density. It is still possi- 
ble to form stars in the field, but cells that are influenced by a 
SN remnant are favoured. Stars which form out of material en- 
riched by a single SN inherit its abundances and therefore show 
an abundance pattern which is characteristic for this particular 
progenitor mass. The most metal-poor stars that form out of 
material which was enriched by only one SN would therefore 



D. Argast et al.: Metal-poor halo stars as tracers of ISM mixing processes 



3 



allow to reconstruct the stellar yields of single core-collapse 
SN, if the average swept up mass and the absolute yield of one 



element were known (Shigeyama & Tsujimoto 1998) 



The SN remnant expansion is the only dynamical process 
taken into account in our model. Therefore, this model has the 
least possible mixing of the ISM. This is the opposite limit as in 
the case of closed box models, which assume a complete mix- 
ing of the ISM at all times and are therefore not able to explain 
the scatter seen at low metallicities. We continue our calcula- 
tion up to an averaged iron abundance of [Fe/H] = —1.0. At 
this metallicity SN events of Type la, which we have not in- 
cluded in our model, would start to influence the ISM. 

3. SN II yields and their uncertainties 

In the present investigation we make use of the nucleosynthesis 
results by Thielemann et al. (1996) and Nomoto et al. (1997). 
Here we want to give a short summary of the key features to- 
gether with an assessment of the uncertainties by comparing 
with available independent calculations. The synthesized ele- 
ments form three different classes which are sensitive to dif- 
ferent aspects of the stellar models and supernovae explosion 
mechanism: (1) stellar evolution, (2) stellar evolution plus the 
explosion energy, and (3) details of the explosion mechanism 
which includes aspects of stellar evolution determining the size 
of the collapsing Fe-core. Due to reaction equilibria obtained in 
explosive burning, the results do not show a stro ng sen sitivity 
to the applied reaction rate library (Hoffman et al. 1999 ). 

1: The abundances of C, O, Ne, and Mg originate from 
the unaltered (essentially only hydrostatically processed) C- 
core and from explosive Ne/C-burning. They are mainly de- 
pendent on the structure and zone sizes of the pre-explosion 
models resulting from stellar evolution. These zones and there- 
fore the amount of ejected mass varies strongly over the pro- 
genitor mass range. O, Ne, and Mg vary by a factor of 10-20 
between a 13 Mq and a 25 Mq progenitor star. This behaviour 
can vary with the treatment of stellar evolution and is strongly 
related to the amount an d meth od of mixing in unstable lay- 
ers. Woosley & Weaver ( 1995 ) employ the Ledoux criterion 
with semiconvection for Schwarzschild-unstable but Ledoux- 



stable layers. Nomoto & Hashimoto ( 1988 ) make use of the 
Schwarzschild criterion for convection (neglecting composi- 
tion gradients) which ensures mixing over more extended re- 
gions than the Ledoux criterion. The Schwarzschild criterion 
causes larger convective cores (see also Chieffi et al. 1998) 
which leads to larger le O, 20 Ne, and 24 Mg yields, the latter 
being also dependent on the 12 C(a, 7) rate (Langer & Henkel 



1995). In addition, it is important to know the mixing veloc- 



ity in unstable regions. Recent calculations by Umeda et al. 



(1999), within the diffusion approximation for mixing (Spruit 



1992; Saio & Nomoto 1998) but with a remaining free parame- 
ter - permitted to vary between and 1 - show that the Woosley 



tional instabilities lead to additional mixing which can bring 
the models making use of the Ledoux criterion closer to those 
evolved with the Schwarzschild criterion and the compositions 
closer to those obtained with instantaneous mixing (high mix- 
ing velocities). Thus, the amount of mixing (being influenced 
by the mixing criterion utilized, the mixing velocity, and ro- 
tation) determines in stellar evolution the size of the C/O core. 
While the yield of O can be fixed with a combination of the still 



uncertain C(a, 7) rate (Buchmann 1996, 1997) and a mixing 



description, the yields of Ne and Mg depend on the extent of 
mixing. Recent galactic che mical (but not dynami c) evo lution 
calculations (Thomas et al. |l998|; Matteucci et al. |l999|; Chi- 



appini et al. 1999 ) prefer apparently a larger extent of mixing 
(caused by either of the effects mentioned above) in order to 
reproduce the observed Mg in low metallicity stars. 

2: The amount of mass for the elements Si, S, Ar and Ca, 
originating from explosive O- and Si-burning, is almo st the 
same for all massive stars in the Thielemann et al. (1996) mod- 



els. They do not show the strong progenitor mass dependence 
of C, O, Ne, and Mg. Si has some contribution from hydrostatic 
burning and varies by a factor of 2-3. Thus, the first set of ele- 
ments (C, O, Ne, Mg) tests the stellar progenitor models, while 
the second set (Si, S, Ar, Ca) tests the progenitor models and 
the explosion energy, because the amount of explosive burn- 
ing depends on the structure of the model plus the energy of 
the shock wave which passes through it. Present models make 
use of an artificially induced shock wave via thermal energy 
deposition (Thielemann et al. 1996) or a piston (Woosley & 



Weaver 1995) with shock energies which lead, after the reduc- 



tion of the gravitational binding of ejected matter, to a given 
kinetic energy. In our models this is an average energy of 10 51 
erg, known from remnant observations, which does not reflect 
possible explosion energy variations as a function of progen- 
itor mass. The apparent underproduction of Ca se en in some 
chemical ev olution calculations (e. g. Tho mas et al. 1998 ; Mat- 
teucci et al. 1999 ; Chiappini et al. 1999 ) could apparently be 
solved by a progenitor mass dependent explosion mechanism 
and energy. 

3: The amount of Fe-group nuclei ejected (which includes 
also one of the so-called alpha elements, i.e. Ti) and their rela- 
tive composition depends directly on the explosion mechanism, 
connected also to the size of the collapsing Fe-core. Observa- 
tional checks of individual supernovae are presently still re- 
quired to test the detailed working of a supernova. The present 
situation is still uncertain and depends on Fe-cores from stel- 
lar evolution, the supranuclear equation of state and maximum 
neutron star mass, related to the total amount of gravitational 
binding energy release of the collapsed protoneutron star, the 
resulting total amount of neutrinos, and the time release (lu- 
minosity), dependent on neutrino transport via numerical treat- 



& Weaver ( 1995) results can be reproduced with a small choice Burrows et al. 1995, 199f; Reddy & Prakash 1997; Burrows 



of this parameter of 0.05. A further effect is due to ro tation. 
Whe n also treating rotation cor rectly (Langer et al. 
et al. |1997| ; Meynet & Maeder [l997| ; Heger et al. 



1997; Talon 



ment, convective transport, and opacities (Burrows 1990; Her 
ant et al. |1994| ; Janka & Mu ller |1995[ |1996j ; Keil & Ja nka 



& Sawyer 1998; Mezzacappa et al. 1998; Messer et al. 1998 



1995 



1999), rota- 



Yamada et al. |1999| ; Pons et al. |1999J). Three types of uncer- 
tainties are inherent in the Fe-group ejecta, related to (i) the 
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total amount of Fe (-group) nuclei ejected and the mass cut 
between neutron star and ejecta, mostly measured by 56 Ni de- 
caying to 56 Fe, (ii) the total explosion energy which influences 
the entropy of the ejecta and with it the degree of alpha-rich 
freeze-out from explosive Si-burning and the abundances of 
radioactive 44 Ti as well as 48 Cr, the latter decaying later to 
48 Ti and being responsible for elemental Ti, and (iii) finally 
the neutron richness or Y c =< Z/A > of the ejecta, dependent 
on stellar structure and the delay time between collapse and ex- 
plosion. Y c influences strongly the ratios of isotopes 57/56 in 
Ni (Co, Fe) and the overall elemental Ni/Fe ratio. The latter be- 
ing dominated by 58 Ni and 56 Fe. The position of the mass cut 
has also a side effect (besides determining the total amount of 
56 Ni/Fe), it influences the ratio of abundances from alpha-rich 
freeze-out and incomplete Si-burning, affecting in this way the 
abundances of the elements Mn ( 55 Co decay), Cr ( 52 F e deca y) 
and Co ( 59 Cu decay) as discussed in Nakamura et al. ( 1999| ). 



There is limited direct observational information from 
individual supernovae with known progenitors (SN 1987 A, 
SN 1993J, 1997D?, 1996N?, 19941) and possible hypernovae 
(SN1997ef, 1998bw), leading to direct O, Ti or Fe (Ni) ob- 
servations (e.g. Shigeyama & Nomoto |l 990^ Iwamoto et al. 



1994, 1998, 1999a, 1999b; Turatto et al. 1998; Sollerman et al 



1998; Kozma & Fransson 1998; Bouchet et al. 1991 ; Suntzeff 



et al. 1992). As explosive nucleosynthesis calculations cannot 
presently rely on self-consistent explosion models, the position 
of the mass cut is in all cases an assumption and has mostly 
been normalized to observations of SN 1987 A. Whether there 
is a decline in Fe-ejecta as a function of progenitor mass (as 
assumed in Thielemann et al. 1996) or actually an increase 



(Woosley & Weaver 1995) or a more complex rise, maximum 



and decline (Nakamura et al. |1999J) is not r eally understood. 
Thus, the results by Thielemann et al. ( 1996 ) utilized here are 
showing the correct IMF integrated behaviour of e.g. Si/Fe, but 
one has to keep in mind that e.g. O/Fe, Mg/Fe, Si/Fe, Ca/Fe 
yields of individual supernovae could be quite uncertain and 
even show an incorrect progenitor mass dependence or a larger 
scatter than (yet unknown) realistic models. Ratios within the 
Fe-group (like e.g. Ni/Fe) have been obtained by mass cut po- 
sitions which reproduce the solar ratios. Thus, the theoreti- 
cal yields might show already the average values and a much 
smaller scatter than some observations (see e.g. Henry 1984). 
Later work attempted to choose mass cuts in order to represent 
some specific element trends like e.g. in Cr/Fe, Co/Fe or Mn/Fe 
(Nakamura et al. |1999(). 



In general we should keep in mind that as long as the ex- 
plosion mechanism is not completely and quantitatively under- 
stood yet, one has to assume a position of the mass cut. Depen- 
dent on that position, which is a function of explosion energy 
and the delay time between collapse and final explosion, the to- 
tal amount of Fe-group matter can vary strongly, Ti-yields can 
vary strongly due to the attained explosion energy and entropy, 
and the ejected mass zones will have a variation in neutron ex- 
cess which automatically changes relative abundances within 
the Fe-group, especially the Ni/Fe element ratio. 



4: r-Process Yields. SNe II have long been expected to 
be the source of r-process elements. Some recent calculations 
seemed to be able to reproduce the solar r-process abundances 
well in the high entropy neutrino wind, emitted from the hot 
protoneutron star after the SN II explosion (Takahashi et al. 
1994; Woosley et al. 1994). If the r-process originates from 



supernovae, a specific progenitor mass dependence has to be 
assumed in order to reproduce the r-process abundances in 
low m etallicity stars as a f unction of [Fe/H] (Mathews et al. 
1992 ; Wheeler et al. 1998 ). Such a "hypothetical" r-process 



yield curve has been constructed by Tsujimoto & Shig eyama 
( 1998 ), in agreement with i deas of Ishimaru & Wanajo ( 1999| ) 
and Travaglio et al. ( 1999 ), and is used in the present galac- 
tic evolution calculation. However, we should keep in mind 
that present-day supernova models have difficulties to repro- 
duce the entropies required for such abundance calculations. In 
addition, they could exhibit the incorrect abu ndance features of 
lighter r-process nuclei (Freiburghaus et al. 1999a), we know 



by now that at least two r-process sources have to contribute to 
the solar r-process abundances (Wasserburg et al. 1996; Cowan 
et al. |1999), and that possible other sources exist (Freiburghaus 



et al. |1999bj). A larger scatter in the r/Fe ratio in low metallicity 
stars than predicted by the constructed supernova yields would 
also indicate the need of such another r-process source. 

4. Observational Data 

As shown in Table |J the observational data were selected 
from various high-resolution studies. Elements marked with 
"c" were adjusted to the same solar abundance scale by S. G. 
Ryan (Ryan et al. 1996 and private communication), entries 
marked with "x" remained unaltered. Unmarked entries were 
not observed in the corresponding study. All observations of 
very metal -poor stars published after 1989 were taken into ac- 
count. In the case where multiple observations of stars exist, we 
used the most recent data. If they were published in the same 
year or after 1995, the values were averaged. 

Typical abundance errors given for these observations are 
about 0.1 dex (see Ryan et al. 1996| ). These may not fully ac- 
count for systematic errors, which could, e.g., be caused by the 
choice of the employed stellar atmosphere models and param- 
eters, or by the employed solar abundance values. If there are 
systematic offsets between different subsamples, this could en- 
hance the scatter in the combined sample. 

Special attention has to be paid to Cr, Mn and O. In the case 
of Cr and Mn, data published after 1995 show a decrease in the 
[Cr/Fe] and [Mn/Fe] ratios for lower metallicities. Such trends 
are not present in older observations. Regardless of this, it is 
not possible to reproduce metallicity-dependent trends with our 
model, since we only use the stellar yields of Thielemann et al. 



(1996 ), which assume constant solar progenitor star metallicity. 



A possible explanation of the [Cr/Fe] and [Mn/Fe] trends is the 
dependency of stellar yields on the metallicity of the progeni- 



tor star of a SN Type II (Samland 1997) or a progenitor mass 



de pende nt mixing of SN II yields with the ISM (Nakamura et 
al.|l999|). 
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Table 1. Reference list of the observational data. Elements marked with "c" were adjusted to the same solar abundance scale by 
S. G. Ryan (Ryan et al. 1996| and private communication), unaltered observations are marked with "x". Unmarked entries were 
not observed in the corresponding study. Newer publications (after 1995) are listed separately. 



Author 



O 



Ma 



Si 



Ca 



Cr 



Mn 



Fc 



Ni 



Eu 



Gratton ( 1989| ) 
Magain ( 198S| ) 
Molaro & Bonifacio ( 
Molaro & Cast elli (|T990| ) 
Peterson et al 
Zhao & Maga in (B99 (ii 
Bessell et al. ( Jl 99T| 
Gratton & Sneden ( 11991a ) 
Gratton & S neden ( 
Ryan et al. (|l99l|) 



1991b) 



Spiesman & Wallers tein ( 199 1 ) 
Spite & Spite ( |l99l|) 
Francois et al. (|1993|) 



Norris et al. (1993) 



Bever idge & Sneden ( |1994| ) 
King (|l994|) 



Nissen et al. (119941) 
Primas et al. Ql994h 
Sneden et al. (119941) 



Fuhrmann et al. ( 1995 ) 
McWilliam et al. (119951) 



B alachandr an & C arney ( 1 996 ) 
Ryan et al. ( |l996|) 
Israelian et al. (|1998|) 



Jehin et al. (1999) 



In the data set of Israelian et al. ( 1998 ), the [O/Fe] ratio 
increases at lower metallicities. This behaviour is apparently 
different from previous determinations and is not yet under- 
stood, but almost identical results were publ ished by a differ- 
ent group of observers (Boesgaard et al. 1999 ). The a-elements 
oxygen and magnesium are produced mainly during the hydro- 
static burning phase of a high mass star and are only slightly 
affected by the actual explosion. Furthermore, their yields de- 
pend in almost the same way on the mass of the progenitor star 
(Thielemann et al. 1996; Woosley & Weaver 1995). Thus ob- 



servations of Mg for the objects discussed in the O abundance 
determinations by Israelian et al. (1998) and Boesgaard et al. 
(1999) would be very interesting. 



5. Results 

5.7. Chemical Mixing of the Halo ISM 

Starting with an ISM of primordial abundances, a first gen- 
eration of ultra metal-poor stars is formed in the model. Af- 
ter the first high-mass stars exploded as SNe of Type II, the 
halo ISM is dominated by local inhomogeneities, since the SN 
events are spatially well separated and no mixing has yet oc- 
curred. In a very metal-poor medium, a single SN event heav- 
ily influences its surroundings, so that its remnant shows the 



element abundance pattern produced by that particular core- 
collapse SN. Stars born out of this enriched material therefore 
inherit the same abundances. Since SNe of different progenitor 
masses have different stellar yields, stars formed out of an in- 
completely mixed ISM show a great diversity in their element 
abundances. 

Fig. [j] shows a cut through the computed volume, giving 
the density distribution of the halo ISM at four different times, 
ranging from the unmixed to the well mixed stage. Each panel 
shows a different mean metallicity and has a lateral length of 
2.5 kpc. The density of the ISM ranges from 6 • 10 -7 Mq/pc 3 
in the inner part of a remnant to 0.04Mq/pc 3 in the densest 
clouds, which is about the gas density of the solar neighbour- 



hood (Binney & Tremaine 1987) 



Upper left: After 82 Myr (see below for the scaling of the 
time units with the assumed SFR), the mean metallicity of the 
halo ISM is [Fe/H] = —3.5. Most of the volume was not yet 
affected by SN events, which can be seen as bright patches in 
the otherwise homogeneously distributed ISM. In the inner part 
of the remnants most of the gas has been swept up and con- 
densed in thin shells, which show up as dark, ring-like struc- 
tures. The regions influenced are well separated and no or only 
slight mixing on a local scale has taken place. The halo ISM 
has to be considered unmixed and is dominated by local inho- 
mogeneities. Stars forming in the neighbourhood of such inho- 
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mogeneities show an abundance pattern which is determined 
by the ejecta of a single SN or a mixture of at most two to three 
SN events. 

Upper right: Beginning of the transition from the unmixed 
to the well mixed halo ISM, after 170 Myr at a mean metallicity 
of [Fe/H] = —3.0. The separation of single SN remnants has 
become smaller and a higher number of remnants may overlap. 
The abundance pattern in overlapping shells still shows a great 
diversity but is closer to the IMF averaged element abundance 
than at the completely unmixed stage. 

Lower left: After 430 Myr, at a mean metallicity of [Fe/H] 
= —2.5, enough SN events have occurred to pollute the mass 
of the whole ISM twice. Nevertheless there are still cells which 
were not influenced by a S, showing therefore primordial abun- 
dances. 

Lower right: After 1430 Myr, a mean metallicity of [Fe/H] 
= —2.0 is reached. At this time, no patches with primordial 
abundances exist and the ISM starts to be chemically well 
mixed. The once homogeneous, primordial medium has now 
a lumpy structure with large density fluctuations and shows 
mainly an IMF averaged abundance pattern. 

Our calculations confirm that the incomplete mixing of the 
ISM during the halo formation plays a significant role in the 
early enrichment of the metal-poor gas. This can be seen in the 
[El/Fe] ratios of the considered elements, shown in Fig. ^. The 
small, filled squares show the [El/Fe] ratio of single model- 
stars. For comparison, observed metal-poor stars are repre- 
sented by open squares. In the case of O, Cr and Mn, obser- 
vations taken before 1995 were marked with open triangles. 
This was done to highlight possible trends in the [El/Fe] ra- 
tio of these elements. If multiple observations of a metal-poor 
star occurred and the abundances had to be averaged, open di- 
amonds were used (see Sect. |]for details). 

The numbered circles show the [El/Fe] and [Fe/H] ratios 
of a single SN II with the indicated mass. The [Fe/H] ratio in 
the shell of the SN remnants are determined by the mass of the 
exploding star and the 5 x 10 4 Mq of gas that are swept up by 
the shock front. This ratio represents the maximal abundance 
a single SN II can produce and determines the position of the 
open circles on the [Fe/H] axis in Fig. ^. If the swept up ma- 
terial subsequently mixes with the surrounding medium, this 
abundance decreases and stars with lower [Fe/H] abundance 
can be formed. On the other hand, the [El/Fe] ratio is deter- 
mined solely by the stellar yields of a SN II and is indepen- 
dent of the mixing mass. For the uncertainties of Fe-ejecta and 
[El/Fe] ratios in individual SNe II see the discussion in Sect. g[ 

We divide the chemical enrichment process in the early 
evolution of the halo ISM into different enrichment stages: 
at metallicities [Fe/H] < —3.0, the ISM is completely un- 
mixed and dominated by local inhomogeneities, originating 
from SN II events. At about [Fe/H] > -2.0 the halo ISM 
shows an IMF averaged abundance pattern and has to be con- 
sidered chemically well mixed. The continuous transition be- 
tween these two phases is marked by incomplete mixing which 



gradually becomes better, as more and more high-mass stars 
explode as SN II. 

The different phases in the enrichment of the halo ISM seen 
in Fig. H can be distinguished in all [El/Fe] - [Fe/H] plots. At 
very low metallicities ([Fe/H] < —3.0), only a few stars exist, 
showing a considerable spread in their [El/Fe] ratios, ranging 
from 0.5 dex in the case of Ni to more than 2 dex in the cases 
of O, Mg and Eu. At this stage, the scatter of the model stars 
is given by the spread in metallicities of the SN models. At the 
end of this early phase, less than 0.5% of the total halo ISM 
mass has been transformed into stars. 

At a metallicity of -3.0 < [Fe/H] < -2.0, the SN rem- 
nants start to overlap and a first, incomplete mixing occurs. 
New stars form out of material which was influenced by sev- 
eral SNe of different masses. Therefore, they do not show the 
typical abundance pattern of a single SN, but show an average 
of the SNe which contributed to the enrichment of the local 
ISM. The spread in the metallicities gradually decreases from 
[Fe/H] = —3.0 to —2.0, reflecting the ongoing mixing process 
as more and more SNe pollute the ISM. At the beginning of the 
well-mixed phase, star-formation has consumed about 2.5% of 
the total halo ISM mass. 

This late phase is characterized by a well mixed ISM and 
begins at [Fe/H] > —2.0, where the abundance scatter in the 
model is reduced to a third of its initial value. At this stage, 
the whole volume considered was influenced several times by 
SN events. This leads to an IMF averaged [El/Fe] ratio in the 
ISM which is the same as predicted by simple 1-zone models 
and which is observed in stars with metallicity [Fe/H] > —2.0. 
Even in an enriched medium, however, a SN event will still 
have an influence on its neighbourhood, although the change in 
the abundance pattern will not be as prominent as in a very 
metal-poor medium. This explains the fact, that even in the 
well-mixed case, the [El/Fe] ratios show a certain dispersion. 
At [Fe/H] = -1.0, about 8% of the total halo ISM mass has 
been used to form stars. 

The onset of SNe of Type la marks the beginning of a third 
phase in the chemical enrichment history of the galaxy, but we 
will not consider this phase any further. Note, that the mean 
[El/Fe] ratios at [Fe/H] = —1.0 in our model do not have to 
be equal to zero (by definition the solar metallicity), since we 
have neglected the influence of SN la events. 

To quantify the enrichment of the ISM we introduce the 
concept of the polluted mass M po \\ in a unit volume at time r. 
It is defined as the total mass which gets polluted by TVsn iso- 
lated SNe, where jVsn is the number of SNe in the unit volume 
that occurred during the elapsed time r. For a constant mixing- 
mass M sw swept up by a SN event, M po ii = Nsn ■ M sw . Since 
the polluted mass is directly proportional to the number of SNe, 
it can become larger than the total ISM mass in the unit vol- 
ume, Mtot- Furthermore, the pollution factor / po n is defined as 
the ratio M po \i/M to t and only depends on iVgN> for fixed M sw 
and M to t- When / po n = 1, enough SNe have contributed to 
the chemical enrichment to theoretically pollute the entire ISM 
in the unit volume, even though there still may be patches of 
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Table 2. Pollution factor and SN Type II frequency. 



[Fe/H] 


/poll 


iVsN [kpc" 3 ] 


r [Myr] 


-3.5 


0.2 


27 


82 


-3.0 


0.7 


83 


170 


-2.5 


2.1 


260 


430 


-2.0 


6.5 


810 


1430 



material which were not yet affected by any SN. A higher pol- 
lution factor results in a better mixing of the halo ISM, decreas- 
ing the local abundance differences and the amount of material 
with primordial abundances. Therefore, the ratio M sw /Mtot 
determines the mixing efficiency in our model, i.e. how many 
SNe in the unit volume are needed to reach a certain value of 
/poll- Given / po n, M BV/ and the mean, IMF integrated iron yield 
(M Fe ) of a typical SN II, the mean metallicity of the ISM then 
is determined by 



[Fe/H] = log 



iVsN • (M Fc ) 
Mtot. 



C = log 



/poll • (-MFc) 



c. 



where C is the solar iron abundance. 

The local inhomogeneities of the ISM begin to disappear 
when most of the gaseous SN remnants start to overlap. This 
is the case when more or less every cloud in the halo was 
influenced at least once by a SN event, i.e. the pollution fac- 
tor is about equal to one. With the adopted mixing mass of 
M sw = 5 x 1O 4 M this is the case at [Fe/H] w -2.8. This 
metallicity gives an upper limit for the end of the early phase 
and the beginning of the transition to the second, well-mixed 
enrichment phase. 

Table ^ shows the pollution factor needed to reach the mean 
metallicities shown in the panels of Fig. |IJ Also shown are the 
corresponding SN frequency Nsn and elapsed time r, which 
depend on our model parameters. Here, the SN frequency is 
defined as the number of SNe per kpc 3 and depends on the 
total ISM mass in the unit volume, whereas the elapsed time 
scales with the average SFR as 



(SFR) 



1.06- 10" 



where 1.06 ■ 10 -4 Mq yr _1 kpc~ 3 is the mean SFR in the unit 
volume in our model. Note that the evolution of the abundance 
ratios as a function of [Fe/H] is independent of the star forma- 
tion timescale and the SFR specified in the model. 

5.2. Comparison with Observations 

The scatter in the [El/Fe] ratios of the model stars as a func- 
tion of [Fe/H] shows the same general trend for every ele- 
ment considered, independent of the individual stellar yields. 
The inhomogeneous mixing of the very metal-poor halo ISM 
at [Fe/H] < —3.0 leads to a scatter in the [El/Fe] ratios of up 
to 1 dex. This scatter continuously decreases for higher metal- 
licities, reflecting the ongoing mixing of the ISM. At [Fe/H] 
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Fig. 3. Normalized nickel-to-iron ratio [Ni/Fe] as function of 
metallicity [Fe/H]. Symbols are the same as in Fig. ||. The scat- 
ter of the model-stars and the halo stars was normalized to unity 
to highlight the enrichment phases of the halo ISM. 



> — 2.0 the model stars show an IMF averaged abundance pat- 
tern with an intrinsic scatter of about 0.1 to 0.2 dex. This be- 
haviour matches the general trend of the observations well, as 
can be seen in Fig. |[ The observations also show a large scatter 
at low metallicities which again decreases for higher [Fe/H], 
with some exceptions, however: the iron-group elements Cr 
and Mn show a strong decrease in the [Cr/Fe] and [Mn/Fe] ratio 
for lower metallicities. This behaviour can not be reproduced 
with our adopted metallicity-independent stellar yields and the 
progenitor-independent mixing mass. 

Compared to the observations, the distribution of [Ni/Fe] 
ratios of the model stars in Fig. ^| shows a scatter that is much 
too small. This is most likely due to the choice of mass cuts in 
the SNe II models, which have been set with the aim to repro- 
duce the average solar [Ni/Fe] ratio. Thielemann et al. ( 1996| ) 
discuss in detail that large variations can easily occur. See also 
the discussion in Sect. ||. We therefore now want to investigate 
whether the sequence of enrichment stages seen in our model 
is similar to the observed evolution of abundance ratios even in 
cases when the employed yields may be incorrect. 

To this end we have normalized the scatter in [Ni/Fe] of the 
model-stars at low [Fe/H] to unity, and have similarly renor- 
malized the range of values for the observed stellar [Ni/Fe] ra- 
tios to one. The mean values of both distributions were left un- 
changed. The resulting renormalized distributions are shown in 
Fig. ||. The remarkably good agreement of both distributions 
after this procedure indicates that the enrichment history of the 
halo ISM implied by the model is consistent with the data, even 
though the employed Ni yields are not. Based on similar com- 
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Fig. 4. Relative frequency of stars normalized to unity in 
[Si/Fe] bins for three different metallicity ranges (see text). The 
solid line shows observational data, the dashed line the model 
stars. The number of included stars is given in the upper left 
corner. 



parisons, we conclude that the abundance ratio data of most 
elements except Mn and Cr are consistent with the predicted 
enrichment history, and the scatter plots in Fig. ^ can thus be 
used to compare the range of the theoretically predicted nucle- 
osynthesis yields with observations. 

To describe the transition from the metal-poor, unmixed to 
the enriched, well mixed ISM more quantitatively, the relative 
frequency of stars at a given [El/Fe] ratio has been analysed 
for the different enrichment phases. In the case of silicon, this 
detailed enrichment history is shown in Fig. [I| The different en- 
richment phases from [Fe/H] < -3.0, -3.0 < [Fe/H] < -2.0 
and [Fe/H] > —2.0 are given in the panels from top to bottom. 
The solid line shows the relative frequency of observed halo 
stars per [Si/Fe] bin for each enrichment phase and the dashed 
line the relative frequency of computed model stars per bin. To 
account for the effect of observational errors on our data, we 
added a random, normally distributed error with standard de- 
viation 0.1 dex in the [El/Fe] and [Fe/H] ratios to the model 
stars. The bin size in the [Si/Fe] ratio is 0.1 dex for observed 
and computed stars, while the position of the histogram for the 
model stars is shifted by 0.01 dex to the left for better visibil- 
ity. The total number of stars included in the plot is given in the 
upper left corner of each panel, where A^bs and A^ mo( j are the 
number of observed stars and of model stars, respectively. 

In the upper panel, the distributions of both the 22 observed 
and the 4226 model stars show a spread in the [Si/Fe] ratio of 
more than one dex. The distribution of the model stars shows 
two wide, protruding wings and a faint peak at [Si/Fe] « 0.2. 
The "right" wing shows a shallow rise from [Si/Fe] w 1.1 to 
the peak. The "left" wing is not as extended and shows a rather 
steep cutoff at [Si/Fe] w —0.3. This asymmetry is due to the 
nucleosynthesis models of core-collapse SNe, which show a 
more or less constant value of [Si/Fe] w —0.3 for progenitor 
masses in the range of 10 — 13 Mq, as can be seen in Fig. |[ 
The distribution of the halo stars peaks at the same location as 
the model stars but extends only down to [Si/Fe] ks —0.1. We 
attribute this to the poor statistic of the data set, since this gap 
is filled in the middle panel. 

The middle panel of Fig. ^| shows the same distribution for 
the intermediate mixing stage of the ISM. The distribution of 
the model stars now has smaller wings, and peaks at [Si/Fe] 
w 0.3. It is still broader than 1 dex, but the majority of the stars 
fall near the IMF averaged [Si/Fe] ratio. The prominent peak 
is caused by the already well-mixed regions, whereas the broad 
distribution shows that the halo ISM is still chemically inhomo- 
geneous. The peak of the observational sample has shifted by 
about 0.2 dex to the right and lies now at [Si/Fe] w 0.4. Com- 
pared to the prediction of the model, the relative frequency of 
the halo stars is too high in the wings of the distribution and 
too low to the left of the peak. 

The lower panel shows the late stage, where the halo ISM 
is well mixed. The broad wings have completely disappeared 
and only the very prominent peak at the IMF averaged value 
remains. The distributions of the 11 observed stars and the 
370000 model stars are in good agreement. At this metallic- 



D. Argast et al.: Metal-poor halo stars as tracers of ISM mixing processes 



9 



0.6 



0.4 



c 

o 



> 
u 

Q 



O 



o 

tn 



0.2 




0.0 

-4.0 -3.5 -3.0 -2.5 -2.0 -1.5 
[Fe/H] 



Fig. 5. Scatter in [Si/Fe] of the model and observed stars. The 
solid line gives the scatter of the model stars, the dashed and 
dotted lines show the scatter of the model stars, folded with 
an error of 0. 1 and 0.2 dex. Filled squares give the standard 
deviation of the observed stars (see text). 



ity no SN of Type la should have polluted the by now well 
mixed ISM and the metal abundance is high enough to restrict 
the impact of single SN II events on the ISM. 

The most prominent feature which characterizes the differ- 
ent enrichment phases, is the intrinsic scatter in the abundances 
of metal-poor stars. This can be seen in Fig. |5| which shows the 
standard deviation of [Si/Fe] as a function of metallicity [Fe/H] 
for the model and the halo stars. The bin size used to compute 
the standard deviation was 0. 1 dex in metallicity. The solid line 
shows the scatter of the unmodified model stars. The influence 
of observational errors on our data was simulated by adding a 
random, normally distributed error with standard deviation 0. 1 
and 0.2 dex in both [Si/Fe] and [Fe/H]. The resulting scatter 
in dependence of metallicity is given by the dashed and dotted 
lines. In the range of —4.0 < [Fe/H] < —3.0 the scatter has 
a more or less constant value of approximately 0.4 dex. It de- 
clines rather steeply in the range —3.0 < [Fe/H] < —2.0 and 
levels off again at metallicities higher than —2.0, depending on 
the assumed observational errors of 0.0, 0. 1 or 0.2 dex. These 
curves show that for errors in this range the scatter at low metal- 
licities is dominated by the intrinsic differences in the element 
abundances of single stars. 

For comparison, the scatter in the [Si/Fe] ratio of observed 
halo stars is represented by filled squares. The observations 
were binned with a bin size of 0.5 dex to compute the stan- 
dard deviations. To estimate the reliability of their scatter in 
[Si/Fe], we built several new data sets by adding a normally 



distributed random error with standard deviation 0. 1 dex to the 
[Si/Fe] ratio and the metallicity of the stars. For each new data 
set, the standard deviation in the different bins was computed. 
The standard deviation for the results from these artificial data 
sets is given in the plot as 1-er error-bars. The scatter of the 
observed abundance ratios shows nicely the features already 
seen in the curves for the model stars. At the first stage of the 
enrichment, it is approximately constant, followed by a steady 
decline in the intermediate mixing phase. At higher metallici- 
ties, the scatter levels off again. 

Since the scatter in [Si/Fe] at [Fe/H] = -4.0 is about 0.4 
dex the observational errors have little influence on the analysis 
at these low metallicities, unless unknown systematic or confu- 
sion errors were large enough to inflate the scatter at [Fe/H] 
= —4.0 to also about 0.4 dex. On the other hand, observational 
errors do dominate the scatter at metallicities [Fe/H] > —2.0, 
when the halo ISM is well mixed and the intrinsic scatter of the 
stars is negligible compared to the observational errors. 

5.3. Individual Elements and Nucleosynthesis 

Oxygen & Magnesium: As expected, the IMF averaged [El/Fe] 
ratios for O and Mg reproduce the mean abundance of the ob- 
served metal-poor stars nicely. The [O/Fe] ratio seems to be 
slightly too low, whereas [Mg/Fe] is slightly too high, but both 
deviations are smaller than 0.1 dex. No trend in the observa- 
tional data of Mg can be seen and a trend i n O on ly becomes 
visible if the observations of Israelian et al. ( 1998 ) are consid- 
ered. 

An important fact is that the scatter in the data, although in- 
creasing at lower metallicities, does not match the large scatter 
of more than two dex predicted by the stellar yields. Since no 
other mixing effects than the overlapping of SN remnants are 
included in our model, the expected scatter is determined by 
the nucleosynthesis yields. If gas flows and the random motion 
of stars in the halo accelerated the chemical mixing, a smaller 
scatter in the model data would be expected. Even then, the fact 
that the observed stars only show [O/Fe] or [Mg/Fe] ratios cor- 
responding to the stellar yields produced by 18 — 50 Mq SNe 
would remain unexplained. 

If we assume a top-heavy IMF which favoured high-mass 
SNe, this problem could be solved. However, the abundance 
pattern of the other elements should then also reflect this, which 
is not the case. This leaves us with two explanations: Either the 
stellar yields of O and Mg or of Fe are incorrect (or both). 

Since the exact location of the mass cut is not known, 
the actual Fe yields are not very well determined (Woosley 



& Weaver 1995; Thielemann et al. 1996) and direct observa 



tional information which links a progenitor to an ejected Fe 
mass is very limited, maybe with the exception of SN 1987 A 
and 1993J. Otherwise only the IMF integrated Fe-yields are 
constrained, but not necessarily their progenitor mass depen- 
dence (declining, rising, or with a maximum, see Nakamura et 
al. 1999 and Sect. ^|). In order to attain a fit to the observational 
data within our evolution model we would need to decrease the 
Fe yields of the 13 or 15 Mq stars by a factor of six. Without 
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adjusting Fe-yields of the more massive stars, and assuming a 
standard Salpeter IMF, this would increase the IMF averaged 
[O/Fe] and [Mg/Fe] ratio by about 0.3 dex (a factor of two) and 
would therefore result in a much too high IMF averaged value. 
Equally, every other abundance ratio would be affected by this 
change. 

On the other hand, the stellar yields of the a-elements O 
and Mg could be too low for the 13 and 15 Mq progenitors. 
The abundances of O and Mg are mainly determined in the 
hydrostatic burning phases and do not depend heavily on the 
explosion mechanism. Changing the stellar yields of O and Mg 
for the 13 and 15 Mq progenitor stars would therefore require 
to adjust the existing stellar models, which suffer from uncer- 
tainties in the theory of convection and the treatment of rota- 
tion. 

Silicon: The IMF averaged [Si/Fe] ratio and the scatter pre- 
dicted by the stellar yields fit the observations perfectly. The 
decrease in the scatter for higher metallicities and therefore the 
different enrichment phases of the halo ISM are clearly visi- 
ble. At [Fe/H] < —3.0, the scatter in the model points reflects 
the scatter predicted by the stellar yields. During the transition 
from the unmixed to the well mixed ISM, in the range —3.0 < 
[Fe/H] < —2.0, the scatter decreases steadily and leads to the 
IMF averaged [Si/Fe] ratio, which is reached at [Fe/H] > —2.0. 

Calcium: The IMF averaged [Ca/Fe] ratio is about 0.2 dex 
lower than the observed mean for metal-poor stars. If the model 
data is shifted by this value to reproduce the mean of the obser- 
vational data, the scatter for very metal-poor stars and the tran- 
sition to the less metal-poor stars fits the data well. Note that, 
contrary to the other a-elements, the stellar yields of Ca are no 
longer approximately proportional to the mass of the progeni- 
tor. This behaviour becomes more pronounced in the cases of 
Cr, Mn and Ni. 

Chromium & Manganese: The iron-peak elements Cr and 
Mn are both produced mainly during explosive silicon burning 
and show an almost identical, complicated dependence on the 
mass of the progenitor. The IMF averaged [El/Fe] ratio repro- 
duces the observations of the less metal-poor stars well. The 
scatter predicted by the stellar yields is only about 0.8 dex and 
is not as large as for the a-elements. 

A notable feature of the observations is the decrease of the 
[Cr/Fe] and [Mn/Fe] ratios for lower metallicities, seen in the 
newer data. If these trends are real, they can not be reproduced 
by metallicity-independent yields unless one assumes a pro- 
genitor mass dependent amount of mixing with the interstellar 
medium (Nakamura et al. 
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The upper limits of the [El/Fe] ratios of Cr and Mn are 
given by the stellar yields of a 30 Mq SN and correspond to 
the highest [Cr/Fe] and [Mn/Fe] ratios seen in metal-poor stars, 
with the exception of the binary CS 22876-032 which shows an 
unusual high [Mn/Fe] ratio of 0.29 dex. Recent high-signal-to- 
noise, high-resolution data result in a lower [Mn/Fe] value for 
this star, placing it near the IMF averaged value (S. G. Ryan, 
private communication). On the other hand, the lower limits of 
the observed [Cr/Fe] and [Mn/Fe] ratios do not correspond at 
all with those given by the stellar yields, which show a ratio 
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Fig. 6. Metallicity [Fe/H] vs. age of single model stars. 



which is too high by up to 0.5 dex in the Thielemann et al. 



(1990 yields. Only a different choice of mass cuts as a func- 



tion of progenitor mass (Nakamura et al. 1999) or as a function 



of metallicity would be able to rectify this. 

Nickel: The stellar yields of the most important iron-peak 
element besides Fe completely fail to reproduce the observa- 
tions. The scatter of the stellar yields is only about 0.5 dex 
compared to about 1 .2 dex seen in the observational data. Com- 
pared to the mean of the observations, the IMF averaged abun- 
dance of the model-stars is about 0.1 dex to high. The small 
scatter in the [Ni/Fe] ratio originates from an almost constant 
ratio of Ni and Fe yields, which means that in the models Ni 
is produced more or less proportional to iron. To reproduce the 
scatter seen in the observations, Ni would have to depend on 
the progenitor mass differently from Fe. This is a shortfall of 
the employed yields, adjusted to reproduce the average solar 
[Ni/Fe] ratio with the choice of their mass cuts. Varying neu- 
tron excess in the yields, however, can change the [Ni/Fe] ratio 
drastically (Thielemann et al. 1996 and Sect. 

Europium: The r-process element Eu reproduces the scatter 
and the mean of the observational data quite well. But com- 
pared to the a-elements O, Mg, Si and Ca its behaviour is very 
different. The highest [Eu/Fe] ratio is produced by low mass 
SN and the lowest ratio by high mass SN, as required when one 
constructs yields under the assumption that the r-process origi- 
nated from SNe II (see Sect. |J). This is exactly the opposite to 
what is found for the a-elements. Therefore, a top-heavy IMF 
would lead to a steadily increasing [Eu/Fe] ratio, since the most 
massive SNe will explode first (cf. Oxygen & Magnesium). 
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5.4. Age-Metallicity Relation 

Common 1-zone chemical evolution models are based on the 
assumption that the system is well mixed at all times. An 
important consequence of this assumption is a monotonically 
increasing metallicity, which leads to a well defined age- 
metallicity relation. Therefore, it is in principle possible to de- 
duce the age of a star if its metallicity is known. This basic 
assumption of 1-zone chemical evolution models was dropped 
in our stochastic approach. Therefore, it is not surprising that 
the well defined age-metallicity relation has to be replaced by 
a statistical relation. In Fig. ^, the metallicity [Fe/H] of model 
stars is plotted against the time of their formation. Model stars 
are represented by small filled squares. For comparison the 
white line visible in the middle of the black strip shows the 
mean age-metallicity relation, corresponding to the relation 
given by a 1-zone model. As can be seen, there is no clear 
age-metallicity relation at any time. Stars which were formed 
in the first 500 million years show a metallicity ranging from 
[Fe/H] < —4.0 up to [Fe/H] > —2.0 and in one extreme case 
up to [Fe/H] « —1.5. On the other hand, stars with metallic- 
ity [Fe/H] — —2.0 could have formed at any time in the first 
3 x 10 9 years. Taking these huge uncertainties into account, it 
is no longer possible to speak of a well-defined age-metallicity 
relation. 

While the mean [Fe/H] abundance increases about linearly 
with time, the scatter in Fig. ^ again reflects the different en- 
richment phases of the ISM. The steep rise at early times marks 
the metal-poor and chemically inhomogeneous stage of the 
halo ISM. Contrary to this first phase, where SN events dom- 
inated the metal-poor ISM locally and the enrichment of iso- 
lated clouds could be very efficient, the late stage is charac- 
terized by a well mixed ISM and therefore an inefficient en- 
richment, which is reflected by the slow increase of the age- 
metallicity "relation" at later times. 

5.5. Ultra Metal-poor Stars 

From our calculations, we can deduce the number of metal- 
poor stars which we expect to observe in different metallicity 
bins. The normalized distribution is shown in Fig. for bin 
sizes of 0.2 dex and 1.0 dex. The number distribution is ap- 
proximately a power law with slope 0.7 for [Fe/H] < —3.5 
and slope 0.9 for [Fe/H] > -3.0. Also plotted are the ob- 
served data from a homogeneous intermediate resolution sam- 
ple of Ryan & Norris (1991) and the combined high resolution 
data from Table |lj both rebinned to the large bin size (see Ta- 
bles U and |] for the numerical values). It is possible that the 
high resolution sample is incomplete in the range [—3.0, —2.0], 
whereas the intermediate resolution sample might have a short- 
age of stars with [Fe/H] < —3.0. Table [?]also lists the num- 
ber of model stars expected in the three 1 .0 dex bins, normal- 
ized such that the number of model stars in the range — 4.0 < 
[Fe/H] < —3.0 is equal to the number of observed halo stars 
in the high-resolution sample in this metallicity range. As can 
be seen, in this case we expect 5 ± 2 model stars with —5.0 < 



Table 3. Top: Relative frequency of stars in the homogeneous 
intermediate resolution survey of Ryan & Norris (1991), the 
combined high resolution data from Table [j] and our model, 
binned with binsize 1 dex. 

Bottom: Absolute numbers. The last row gives the number of 
stars per bin which we expect to be present, if our model gives 
a fair representation of the halo metallicity distribution. The 
number of model-stars is normalized to the number of stars 
in the range —4.0 < [Fe/H] < —3.0 in the high resolution 
sample. No star was detected with confirmed [Fe/H] < —4.0, 
in contrast to the 5 ± 2 stars predicted by the model. 



[Fe/H] 


-3.0,-2.0] 


[-4.0, -3.0] [- 


5.0, -4.0] 


Ryan & Norris 


0.943 


0.057 


0.000 


High resolution 


0.654 


0.346 


0.000 


Model 


0.865 


0.115 


0.020 


Ryan & Norris 


99 


6 





High resolution 


53 


28 





Expected 


211 ± 15 


28 ±5 


5±2 


Table 4. Relative frequency of model stars, binned with binsize 


0.2 dex. 








[Fe/H] 


N/N tot 


\og{N/N to t) 



[-5.0, -4.8] 


1.92 


1Q -03 


-2.72 


[-4.8, -4.6] 


2.37 


10 -oa 


-2.63 


[-4.6, -4.4] 


3.43 


1Q -03 


-2.46 


[-4.4, -4.2] 


5.10 


1Q -03 


-2.29 


[-4.2, -4.0] 


7.51 


1Q -03 


-2.12 


[-4.0,-3.8] 


1.00 


1Q -02 


-2.00 


[-3.8,-3.6] 


1.31 


1Q -02 


-1.88 


[-3.6,-3.4] 


1.93 


1Q -02 


-1.71 


[-3.4, -3.2] 


2.82 


1Q -Q2 


-1.55 


[-3.2,-3.0] 


4.41 


1Q -Q2 


-1.36 


[-3.0, -2.8] 


6.71 


1Q -Q2 


-1.17 


[-2.8,-2.6] 


1.02 


10 -oi 


-0.99 


[-2.6,-2.4] 


1.47 


10 -oi 


-0.83 


[-2.4, -2.2] 


2.11 


10 -oi 


-0.68 


[-2.2, -2.0] 


3.38 


1Q -oi 


-0.47 



[Fe/H] < —4.0 while the high resolution sample contains none. 
If the ratio of stars in these two metallicity bins for this admit- 
tedly inhomogeneous sample is representative for the Galactic 
halo stars, this would suggest a genuine shortage of the most 
metal-poor stars. In this case a possible solution could be the 
pre-enrichment of the halo ISM by population III stars. It is 
conceivable that these already produced an iron abundance of 
[Fe/H] « —4.0 before they disappeared, leaving only a pre- 
enriched ISM. 



6. Conclusions & Discussion 

We have developed a stochastic model of the early chemical 
enrichment of the halo ISM. The aim of the model is to un- 
derstand the scatter in the [El/Fe] ratios of observed stars at 
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[Fe/H] 



Fig. 7. Metallicity distribution in the model at the end of the 
calculation ([Fe/H]= —1). The number of stars with metallicity 
[Fe/H] is shown (a) in bins of 0.2 dex (b) in bins of 1.0 dex. In 
both cases the total number of stars is normalized to one. The 
solid squares show the intermediate resolution data of Ryan 
& Norris (1991), the solid triangles the combined high resolu- 
tion data from Table |lj both binned with bin size 1.0 dex and 
similarly normalized. There are no observed stars with [Fe/H] 
< -4.0. See Tables |, | and text. 



very low metallicities and the transition to the smaller scatter 
seen at higher metallicities. We computed the evolution of the 
abundances of the a-elements O, Mg, Si and Ca, the iron-peak 
elements Cr, Mn, Fe and Ni and the r-process element Eu and 
investigated the mixing of the halo ISM. 

We divide the enrichment history of the halo ISM into dif- 
ferent phases. At metallicities [Fe/H] < —3.0, the ISM is not 
mixed and is dominated by local abundance inhomogeneities, 
which are caused by individual Type II SNe of different pro- 
genitor masses. The second phase at [Fe/H] > —2.0 is defined 
by a well mixed ISM which shows IMF averaged abundances 
and is too metal-rich to be dominated by single SN events. A 
continuous transition from the first to the second phase occurs 
between -3.0 < [Fe/H] < -2.0. The onset of SN la events 
marks the beginning of a third phase in the enrichment of the 
ISM. 

The different enrichment phases of the halo ISM can be 
distinguished in every considered element, after normalising to 
the stellar yields. On the other hand, Fig. |2| shows that some 
stellar yields reproduce the IMF averaged element-to-iron ratio 
[El/Fe] well, but fail to reproduce the abundance scatter ob- 
served in metal-poor halo stars. Especially the stellar yields of 
the a-elements O and Mg predict stars with a low element-to- 



iron ratio ([O/Fe] or [Mg/Fe] < 0.0), which are not observed. 
Typically, metal-poor halo stars show an overabundance of a- 
elements of about [a/Fe] w 0.3 to 0.4, as can be seen in Fig. |[ 
This is especially troublesome, since an attempt to solve this 
problem would either require a change in the iron yields of 
the 13 and 15 Mg models of up to a factor six, which would 
raise every other mean element abundance by about 0.3 dex (a 
factor of two), or a change in the stellar yields of oxygen and 
magnesium, which are produced mainly during the hydrostatic 
burning phases. 

Recent observations of metal-poor stars show a decrease of 
the [Cr/Fe] and [Mn/Fe] ratios for lower metallicities. It is not 
possible to explain these trends with the metallicity indepen- 
dent stellar yields we have used. Also, the stellar yields of the 
iron-peak element Ni predict a scatter in [Ni/Fe] which is much 
too small compared to the observational data. These problems 
argue strongly for a revision of the theoretical nucleosynthesis 
models and their extension to lower metallicities. 

The unfortunate situation is, however, that there exists no 
theoretical foundation to do so for Fe-group yields as long as 
the supernova explosion mechanism is not understood. Thiele- 
mann et al. ( |1996| ) discussed in detail uncertainties of Fe-group 
yields due to the choice of mass cut, explosion energy and en- 
tropy, as well as the delay time between collapse and explosion, 
affecting also the neutron-richness of matter. Multidimensional 
aspects might add further degrees of freedom (Nagataki et al. 



1997, 1998). Thus, there is an understanding of dependences, 
coincidences of abundance features etc., but at present only ob- 
servational information combined with galactic evolution mod- 
elling like in the present paper or, e.g., Tsujimoto & Shigeyama 



(1998) and Nakamura et al. ( 1999) can try to provide sufficient 



constraints for a further understanding of supernova nucleosyn- 
thesis. 

The advancing enrichment process of the halo ISM can be 
characterized by the pollution factor / po ii, defined as the ratio 
of the mass M po \\ = A*sn ■ Msw polluted by ]Vsn preceding 
SNe and the ISM mass M to t, all in a unit volume. The en- 
richment history of the halo ISM now is mainly determined 
by the mixing efficiency which in turn is fixed by the ratio of 
the mass M sw of swept-up material in a SN event and M tot . 
The more mass a SN sweeps up, the less SN events are needed 
to reach a certain value of / po ii, making the mixing more ef- 
ficient. M sw also determines the average metallicity [Fe/H] of 
the halo ISM for a given pollution factor. A larger swept-up 
mass leads to a lower mean ISM metallicity and vice versa. 
Therefore, the metallicity where the transition from one en- 
richment stage to the next occurs, depends only on the mixing 
efficiency M sw /M to t and not on the SFR. The SFR is only im- 
portant if one is interested in the elapsed time that is needed to 
reach a certain mean metallicity of the halo ISM. 

Figures ^| and ^ support our adopted value of 5 x 10 4 Mq 
of swept up material. If the swept up mass is higher, the mixing 
would be more efficient, resulting in an IMF averaged chemical 
abundance pattern at lower metallicities. This would produce a 
narrower peak in Fig. Q and a steeper slope in Fig. More- 
over, the whole curve in Fig. |] would be shifted towards lower 
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metallicities. On the other hand, a smaller mixing mass would 
reduce the efficiency of the enrichment of the ISM, which could 
be seen in a broader distribution at higher metallicities in Fig. ^| 
and a shallower slope in Fig. |[ 

Standard 1-zone chemical evolution models predict a well 
defined age-metallicity relation, based on the assumption that 
the ISM is well mixed at all times. In our stochastic model 
the chemical inhomogeneity of the halo ISM and therefore 
the scatter in metallicity at any time is much too high to rea- 
sonably establish such a relation. Nevertheless, the steep rise 
with large scatter seen in Fig. || at very early times marks the 
chemically inhomogeneous enrichment phase, while the slow 
increase later-on reflects the well mixed, metal-rich ISM. 

The results of our model also quantify the problem of the 
missing ultra metal-poor stars. From it we have deduced the ex- 
pected number of ultra metal-poor stars (with [Fe/H] < —4.0) 
which should have been observed, normalized to the number of 
halo stars in the combined high-resolution sample with metal- 
licities in the range —4.0 < [Fe/H] < —3.0. We expect about 
5±2 ultra metal-poor stars whereas none was found to date. It is 
possible that Population III stars have caused a pre-enrichment 
of the ISM to [Fe/H] w —4.0 and already have disappeared 
before the onset of the formation of the Galaxy. 

We would like to thank the referee S. G. Ryan for helpful 
and motivating discussions and for providing us with some of 
the observational data. This work was supported by the Swiss 
Nationalfonds. 
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Fig. 1. Cut through the computed volume, showing the density distribution of the halo ISM during the transition from the unmixed 
to the well mixed stage (see text for details). 



16 



D. Argast et al.: Metal-poor halo stars as tracers of ISM mixing processes 




■ ■©■■ © 

H ' 1 ' 1 ' 1 ' h 





o 




□ D 

1 > : ^nK>T?^s^rf M 



1.5 



1.0 



„ 0.5 



" 0.0 



-0.5 



-1.0 



1.5 



1.0 



„ 0.5 



-0.5 



-1.0 



1.5 



1.0 



0.5 



— 0.0 



-0.5 



-1.0 



1.5 



1.0 



,-, 0-5 



— 0.0 



-0.5 



-1.0 



o 




.© 




' D -' $ '.'.■^.'3 ;: ^i 5 ^ 



■©.-:■ i©" 




0. 




-4.5 -4.0 -3.5 -3.0 -2.5 -2.0 -1.5 -4.5 -4.0 -3.5 -3.0 -2.5 -2.0 -1.5 

[Fe/H] [Fe/H] 

Fig. 2. Element-to-iron ratio [El/Fe] of O, Mg, Si, Ca, Cr, Mn, Ni and Eu. Open circles depict [El/Fe] ratios of SN II models of 
the given progenitor mass. Small filled squares represent model stars, open symbols show observed stars. 



